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2D ARRAYS OF JOSEPHSON NANOCONTACTS 
AND NANOGRANULAR SUPERCONDUCTORS 

ABSTRACT 

By introducing a realistic model of nanogranular superconductors (NGS) based on 2D arrays of 
Josephson nanocontacts (created by a network of twin-boundary dislocations with strain fields 
acting as insulating barriers between hole-rich domains), in this Chapter we present some novel 
phenomena related to mechanical, magnetic, electric and transport properties of NGS in un- 
derdoped single crystals. In particular, we consider chemically induced magnetoelectric effects 
and flux driven temperature oscillations of thermal expansion coefficient. We also predict a 
giant enhancement of the nonlinear thermal conductivity of NGS reaching up to 500% when 
the intrinsically induced chemoelectric field (created by the gradient of the chemical potential 
due to segregation of hole producing oxygen vacancies) closely matches the externally pro- 
duced thermoelectric field. The estimates of the model parameters suggest quite an optimistic 
possibility to experimentally realize these promising and important for applications effects in 
non-stoichiometric NGS and artificially prepared arrays of Josephson nanocontacts. 

1. INTRODUCTION 

Inspired by new possibilities offered by the cutting-edge nanotechnologies, the experi- 
mental and theoretical physics of increasingly sophisticated mesoscopic quantum devices, 
heavily based on Josephson junctions (JJs) and their arrays (JJAs), is becoming one of 
the most exciting and rapidly growing areas of modern science (for reviews on charge and 
spin effects in mesoscopic 2D JJs and quantum-state engineering with Josephson devices, 
see, e.g., Newrock et al. 2000, Makhlin et al. 2001, Krive et al. 2004, Sergeenkov 2006, 
Beloborodov et al. 2007). In particular, a remarkable increase of the measurements tech- 
nique resolution made it possible to experimentally detect such interesting phenomena as 
flux avalanches (Altshuler and Johansen 2004) and geometric quantization (Sergeenkov 
and Araujo-Moreira 2004) as well as flux dominated behavior of heat capacity (Bourgeois 
et al. 2005) both in JJs and JJAs. 

Recently, it was realized that JJAs can be also used as quantum channels to transfer 
quantum information between distant sites (Makhlin et al. 2001, Wendin and Shumeiko 
2007) through the implementation of the so-called superconducting qubits which take 
advantage of both charge and phase degrees of freedom. 

Both granular superconductors and artiflcially prepared JJAs proved useful in studying 
the numerous quantum (charging) effects in these interesting systems, including Coulomb 
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blockade of Cooper pair tunneling (lansity et al. 1988), Bloch oscillations (Haviland et al. 
1991), propagation of quantum ballistic vortices (van der Zant 1996), spin-tunneling re- 
lated effects using specially designed 5'F5'-type junctions (Ryazanov et al. 2001, Golubov 
et al. 2002), novel Coulomb effects in SINIS-tjpc nanoscale junctions (Ostrovsky and 
Feigel'man 2004), and dynamical AC reentrance (Araujo-Moreira et al. 1997, Barbara et 
al. 1999, Araujo-Moreira et al. 2005). 

At the same time, given a rather specific magnetostrictive (Sergeenkov and Ausloos 
1993) and piezomagnetic (Sergeenkov 1998b, Sergeenkov 1999) response of Josephson sys- 
tems, one can expect some nontrivial behavior of the thermal expansion (TE) coefficient in 
JJs as well (Sergeenkov et al. 2007). Of special interest are the properties of TE in applied 
magnetic field. For example, some superconductors like Bai-xK^BiO^, BaPhxBii^xO^ 
and La2-xSrxCu04 were found (Anshukova et al. 2000) to exhibit anomalous tempera- 
ture behavior of both magnetostriction and TE which were attributed to the field-induced 
suppression of the superstructural ordering in the oxygen sublattices of these systems. 

The imaging of the granular structure in underdoped Bi2Sr2CaCu20s+s crystals 
(Lang et al. 2002) revealed an apparent segregation of its electronic structure into super- 
conducting domains (of the order of a few nanometers) located in an electronically distinct 
background. In particular, it was found that at low levels of hole doping (5 < 0.2), the 
holes become concentrated at certain hole-rich domains. (In this regard, it is interesting 
to mention a somewhat similar phenomenon of " chemical localization" that takes place in 
materials, composed of atoms of only metallic elements, exhibiting metal-insulator tran- 
sitions, see, e.g., Gantmakher 2002.) Tunneling between such domains leads to intrinsic 
nanogranular superconductivity (NGS) in high-Tc superconductors (HTS). Probably one 
of the first examples of NGS was observed in Y Ba2Cu^Oi-5 single crystals in the form of 
the so-called "fishtail" anomaly of magnetization (Daeumling et al. 1990). The granular 
behavior has been related to the 2D clusters of oxygen defects forming twin boundaries 
(TBs) or dislocation walls within CuO plane that restrict supercurrent flow and allow 
excess flux to enter the crystal. Indeed, there are serious arguments to consider the TB 
in HTS as insulating regions of the Josephson SIS-type structure. An average distance 
between boundaries is essentially less than the grain size. In particular, the networks of 
localized grain boundary dislocations with the spacing ranged from lOnm to lOOnm have 
been observed (Daeumling et al. 1990) which produce effectively continuous normal or 
insulating barriers at the grain boundaries. It was also verified that the processes of the 
oxygen ordering in HTS leads to the continuous change of the lattice period along TB 
with the change of the oxygen content. Besides, a destruction of bulk superconductivity 
in these non-stoichiometric materials with increasing the oxygen deficiency parameter 5 
was found to follow a classical percolation theory (Gantmakher et al. 1990). 

In addition to their importance for understanding the underlying microscopic mecha- 
nisms governing HTS materials, the above experiments can provide rather versatile tools 
for designing chemically-controlled atomic scale JJs and JJAs with pre-selected properties 
needed for manufacturing the modern quantum devices (Sergeenkov 2001, Araujo-Moreira 
et al. 2002, Sergeenkov 2003, Sergeenkov 2006). Moreover, as we shall see below, NGS 
based phenomena can shed some light on the origin and evolution of the so-called para- 
magnetic Meissner effect (PME) which manifests itself both in high-Tc and conventional 
superconductors (Geim et al. 1998, De Leo and Rotoli 2002, Li 2003) and is usually 
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associated with the presence of vr-junctions and/or unconventional ((i-wave) pairing sym- 
metry. 

In this Chapter we present numerous novel phenomena related to the magnetic, elec- 
tric, elastic and transport properties of Josephson nanocontacts and NGS. The paper is 
organized as follows. In Section 1, a realistic model of NGS is introduced which is based 
on 2D JJAs created by a regular network of twin-boundary dislocations with strain fields 
acting as an insulating barrier between hole-rich domains (like in underdoped crystals). 
In Section 2, we consider some phase-related phenomena expected to occur in NGS, such 
as Josephson chemomagnetism and magnctoconcentration effect. Section 3 is devoted to 
a thorough discussion of charge-related polarization phenomena in NGS, including such 
topics as chemomagnetoelectricity, magnetocapacitance, charge analog of the "fishtail" 
(magnetization) anomaly, and field-tuned weakening of the chemically-induced Coulomb 
blockade. In Section 4 we present our latest results on the infiuence of an intrinsic chem- 
ical pressure (created by the gradient of the chemical potential due to segregation of hole 
producing oxygen vacancies) on temperature behavior of the nonlinear thermal conduc- 
tivity (NLTC) of NGS. In particular, our theoretical analysis (based on the inductive 
model of 2D JJAs) predicts a giant enhancement of NLTC reaching up to 500% when 
the intrinsically induced chemoelectric field = ^^A* closely matches thermoelectric 
field Et — 5'tVT. And finally, by introducing a concept of thermal expansion (TE) of 
Josephson contact (as an elastic response of J J to an effective stress field), in Section 5 we 
consider the temperature and magnetic field dependence of the TE coefficient a{T, H) in 
a small single J J and in a single plaquette (a prototype of the simplest JJA). In particular, 
we found that in addition to expected field oscillations due to Praunhofer-like dependence 
of the critical current, a of a small single junction also exhibits strong fiux driven tem- 
perature oscillations near Tc- The condition under which all the effects predicted here 
can be experimentally realized in artificially prepared JJAs and NGS are also discussed. 
Some important conclusions of the present study are drawn in Section 6. 

2. MODEL OF NANOSCOPIC JOSEPHSON JUNCTION ARRAYS 

As is well-known, the presence of a homogeneous chemical potential n through a single 
J J leads to the AC Josephson effect with time dependent phase difference d(f)/dt — ji/h. 
In this Section, we will consider some effects in dislocation induced JJ caused by a local 
variation of excess hole concentration c(x) under the chemical pressure (described by 
inhomogeneous chemical potential /Li(x)) equivalent to presence of the strain field of 2D 
dislocation array e(x) forming this Josephson contact. 

To understand how NGS manifests itself in non-stoichiometric crystals, let us invoke 
an analogy with the previously discussed dislocation models of twinning-induced super- 
conductivity (Khaikin and Khlyustikov 1981) and grain-boundary Josephson junctions 
(Sergeenkov 1999). Recall that under plastic deformation, grain boundaries (GBs) (which 
are the natural sources of weak links in HTS), move rather rapidly via the movement of 
the grain boundary dislocations (GBDs) comprising these GBs. At the same time, ob- 
served (Daeumling et al. 1990, Lang et al. 2002, Yang et al. 1993, Moeckley et al. 
1993) in HTS single crystals regular 2D dislocation networks of oxygen depleted regions 
(generated by the dissociation of < 110 > twinning dislocations) with the size do of a 
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few Burgers vectors, forming a triangular lattice with a spacing d > do ranging from 
lOnm to lOOnm, can provide quite a realistic possibility for the existence of 2D Josephson 
network within CuO plane. Recall furthermore that in a d-wave orthorhombic YBCO 
crystal TBs are represented by tetragonal regions (in which all dislocations are equally 
spaced by d^ and have the same Burgers vector a parallel to |/-axis within CuO plane) 
which produce screened strain fields (Gurevich and Pashitskii 1997) e(x) = eoe"''''/'^" with 
I X I = >/a;2 + y'^. 

Though in YBa2Cu^O-j_s the ordinary oxygen diffusion D — DQt"^''-/^'^'^ is extremely 
slow even near Tc (due to a rather high value of the activation energy Ud in these materials, 
typically Ud — leV), in underdoped crystals (with oxygen-induced dislocations) there is a 
real possibility to facilitate oxygen transport via the so-called osmotic (pumping) mecha- 
nism (Girifalco 1973, Sergeenkov 1995) which relates a local value of the chemical potential 
(chemical pressure) //(x) = //(O) -|- • x with a local concentration of point defects as 
follows c(x) = e~^'^^^^^^'^ . Indeed, when in such a crystal there exists a nonequilibrium 
concentration of vacancies, dislocation is moved for atomic distance a by adding excess 
vacancies to the extraplane edge. The produced work is simply equal to the chemical 
potential of added vacancies. What is important, this mechanism allows us to explicitly 
incorporate the oxygen deficiency parameter 6 into our model by relating it to the excess 
oxygen concentration of vacancies Cv = c(0) as follows 5 = 1 — As a result, the chemical 
potential of the single vacancy reads /i„ = /i(0) = —ksT \og{l — 5) ~ ksTd. Remarkably, 
the same osmotic mechanism was used by Gurevich and Pashitskii (1997) to discuss the 
modification of oxygen vacancies concentration in the presence of the TB strain field. In 
particular, they argue that the change of e(x) under an applied or chemically induced 
pressure results in a significant oxygen redistribution producing a highly inhomogcncous 
filamentary structure of oxygen-deficient nonsuperconducting regions along GB (Moeck- 
ley et al. 1993) (for underdoped superconductors, the vacancies tend to concentrate in the 
regions of compressed material). Hence, assuming the following connection between the 
variation of mechanical and chemical properties of planar defects, namely /i(x) = i^Qoe(x) 
(where is an effective atomic volume of the vacancy and K is the bulk elastic mod- 
ulus), we can study the properties of TB induced JJs under intrinsic chemical pressure 
V/i (created by the variation of the oxygen doping parameter 5). More specifically, a 
single SIS type junction (comprising a Josephson network) is formed around TB due 
to a local depression of the superconducting order parameter A(x) oc e(x) over distance 
do producing thus a weak link with (oxygen deficiency 5 dependent) Josephson coupling 
J{5) = e(x)Jo = Jo(5)e"l''l/*' where Jo{S) = cqJo = {fJ'v/K^o)Jo (here Jq oc Ao/-R„ with 
Rn being a resistance of the junction). Thus, the present model indeed describes chemi- 
cally induced NGS in underdoped systems (with S 0) because, in accordance with the 
observations, for stoichiometric situation (when 5 ~ 0), the Josephson couphng J (5) ~ 
and the system loses its explicitly granular signature. 

To adequately describe chemomagnctic properties of an intrinsically granular super- 
conductor, we employ a model of 2D overdamped Josephson junction array which is based 
on the well known Hamiltonian 

N N 

n = Y.Mi-cosct>,,) + Y,fL (1) 
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and introduces a short-range interaction between junctions (which are formed around 
oxygen-rich superconducting areas with phases arranged in a two-dimensional (2D) 

lattice with coordinates Xi = {xi,yi). The areas are separated by oxygen-poor insulating 
boundaries (created by TB strain fields e(xjj)) producing a short-range Josephson coupling 
■J'ij = Jo{S)e~^^^^^^'^. Thus, typically for granular superconductors, the Josephson energy 
of the array varies exponentially with the distance Xjj = Xj — Xj between neighboring 
junctions (with d being an average junction size). As usual, the second term in the rhs of 
Eq.(l) accounts for Coulomb effects where = — 2enj is the junction charge with rii being 
the pair number operator. Naturally, the same strain fields e(xij) will be responsible for 
dielectric properties of oxygen-depleted regions as well via the (^-dependent capacitance 
tensor Cij{5) = C[e{xij)]. 

If, in addition to the chemical pressure V/i(x) = KilQVe{x.), the network of supercon- 
ducting grains is under the influence of an applied frustrating magnetic fleld B, the total 
phase difference through the contact reads 

<l>vit) = 4>l + ^(x.,- A n,,) • B + (2) 

where 0°^- is the initial phase difference (see below), n^j — Xjj/| Xjj | with Xj^ = (xj -|- 
Xj)/2, and w = 2Xl{T)+1 with A^^ being the London penetration depth of superconducting 
area and / an insulator thickness which, within the discussed here scenario, is simply equal 
to the TB thickness (Sergeenkov 1995). 

To neglect the influence of the self-fleld effects in a real material, the corresponding 
Josephson penetration length Aj = •>/ $o/2vr/xoJcW must be larger than the junction size 
d. Here is the critical current density of superconducting (hole-rich) area. As we shall 
see below, this condition is rather well satisfied for HTS single crystals. 

Within our scenario, the sheet magnetization M of 2D granular superconductor is 
defined via the average Josephson energy of the array 

r dt f d^x^ . 
<H>^ — 7i(x,i) (3) 



as follows 

M(B.^).-^^, (4) 

where s = 27rd'^ is properly defined normalization characteristic Josephson 

time, and we made a usual substitution j^^^^j Aij{t) ^ J d'^xA{-x,t) valid in the long- 
wavelength approximation (Sergeenkov 2002). 

To capture the very essence of the superconducting analog of the chemomagnetic effect, 
in what follows we assume for simplicity that a stoichiometric sample (with 5 ~ 0) does 
not possess any spontaneous magnetization at zero magnetic field (that is M(0, 0) = 0) 
and that its Meissner response to a small applied field B is purely diamagnetic (that is 
M(B,0) ~ — B). According to Eq.(4), this condition implies = 27rm for the initial 
phase difference with m = 0,±l,±2,... 

Taking the applied magnetic field along the c-axis (and normal to the CuO plane), we 
obtain finally 
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for the chemically-induced sheet magnetization of the 2D Josephson network. Here 
Mo{6) = MS) /Bo with Jo{S) defined earlier, b = B/Bq, and h^, = B^/Bq ~ {kBTT/h)6 
where B^((5) = (/it,r/^)Bo is the chemically-induced contribution (which disappears in 
optimally doped systems with 5 ~ 0), and Bq = ^o/wd is a characteristic Josephson field. 




B/B„ 



Fig. 1: The susceptibility func- 
tion of applied magnetic field for dif- 
ferent values of oxygen deficiency pa- 
rameter: 5 ~ (solid line), 6 = 0.05 
(dashed line), and 6 = 0.1 (dotted 
line) . 




Fig. 2: The oxygen deficiency in- 
duced susceptibility for different 
values of applied magnetic field 
(chemomagnetism) . 



FigHlshows changes of the initial (stoichiometric) diamagnetic susceptibility S) = 
9M(B, 6) 1 9B (solid line) with oxygen deficiency b. As is seen, even relatively small values 
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of 6 parameter render a low field Meissner phase strongly paramagnetic (dotted and dashed 
lines). Figl2] presents concentration (deficiency) induced susceptibility x(B,5)/xo(0) for 
different values of applied magnetic field b = B/Bq including a true chemomagnetic effect 
(solid line). According to Eq.(5), the initially diamagnetic Meissner effect turns paramag- 
netic as soon as the chemomagnetic contribution B^((5) exceeds an applied magnetic field 
B. To see whether this can actually happen in a real material, let us estimate a magnitude 
of the chemomagnetic field B^. Typically (Daeumling et al. 1990, Gurevich and Pashit- 
skii 1997), for HTS single crystals Al(0) ~ 150nm and d ~ lOram, leading to Bo — 0.5T. 
Using T c:^ h/fiy and jc = 10^^ A/m'^ as a pertinent characteristic time and the typical 
value of the critical current density, respectively, we arrive at the following estimate of 
the chemomagnetic field B^((5) ~ 0.5Bq for S = 0.05. Thus, the predicted chemically 
induced PME should be observable for applied magnetic fields B ~ 0.5i?o — 0.25T which 
are actually much higher than the fields needed to observe the previously discussed piezo- 
magnetism and stress induced PME in high-Tc ceramics (Sergeenkov 1999). Notice that 
for the above set of parameters, the Josephson length Aj ~ l^m, which means that the 
assumed here small-junction approximation (with d <^ Aj) is valid and the so-called "self- 
field" effects can be safely neglected. So far, we neglected a possible field dependence of 




Fig. 3: Magnetic field dependence 



0.01 0.1 1 10 of the oxygen vacancy concentration 

^^^0 (magnetoconcentration effect). 

the chemical potential fi^ of oxygen vacancies. However, in high enough applied magnetic 
fields B, the field- induced change of the chemical potential A/it,(B) = /ii,(B) — /it,(0) 
becomes tangible and should be taken into account. As is well-known (Abrikosov 1988, 
Sergeenkov and Ausloos 1999), in a superconducting state A/i^(B) = — M(B)B/?7,, where 
M(B) is the corresponding magnetization, and n is the relevant carriers number density. 
At the same time, within our scenario, the chemical potential of a single oxygen vacancy 
fiv depends on the concentration of oxygen vacancies (through deficiency parameter 6). 
As a result, two different effects are possible related respectively to magnetic field depen- 
dence of /i„(B) and to its dependence on magnetization /i^(M). The former is nothing 
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else but a superconducting analog of the so-called magnetoconcentration effect which was 
predicted and observed in inhomogeneously doped semiconductors (Akopyan et al.l990) 
with field-induced creation of oxygen vacancies Cv(B) = Cv(0) exp^—Afi^CB) / ksT) , while 
the latter results in a "fishtail" -like behavior of the magnetization. Let us start with the 
magnetoconcentration effect. FigJ3] depicts the predicted field-induced creation of oxygen 
vacancies Cv(B) using the above-obtained magnetization M(B,(5) (see Fig.l and Eq.(5)). 
We also assumed, for simplicity, a complete stoichiometry of the system in a zero mag- 
netic field (with Cv(0) = 1). Notice that Cv(B) exhibits a maximum at — 0.23 for 
applied fields B = Bq (in agreement with the classical percolative behavior observed in 
non-stoichiometric Y Ba2CuiOT-s samples (Daeumling et al. 1990, Gantmakher et al. 
1990, Moeckley et al. 1993). Finally, let us show that in underdoped crystals the above- 
discussed osmotic mechanism of oxygen transport is indeed much more effective than a 
traditional diffusion. Using typical YBCO parameters (Gurevich and Pashitskii 1997), 
eo = 0.01, VLq = Oq with = 0.2nm, and K = 115GPa, we have /it,(0) = eoKVlQ ~ ImeV 
for a zero-field value of the chemical potential in HTS crystals, which leads to creation of 
excess vacancies with concentration Cv(0) = e~^"^^'^^^^"^ ~ 0.75 (equivalent to a deficiency 
value of 5(0) ~ 0.25) at T = Tc, while the probability of oxygen diffusion in these materials 
(governed by a rather high activation energy Ud — leV) is extremely low under the same 
conditions because D oc e~^'^/^^'^^ <^ 1. On the other hand, the change of the chemical 
potential in applied magnetic field can reach as much as (Sergeenkov and Ausloos 1999) 
A/i^(B) ~ 0.5meV^ for B = 0.5T, which is quite comparable with the above-mentioned 
zero-field value of /U^(0). Let us turn now to the second effect related to the magnetiza- 



5(0)=0 

5(0)=0.1 

8(0)=0.2 




Fig. 4: A "fishtail" -like behavior of 
magnetization in applied magnetic 
field in the presence of magnetocon- 
centration effect (with field-induced 
oxygen vacancies Cv(B), see Fig. 3) 
for three values of field-free defi- 
ciency parameter: 5(0) ~ (solid 



2 3 4 line), 5(0) = 0.1 (dashed line), and 
B/B^ 5(0) = 0.2 (dotted line). 



tion dependence of the chemical potential /i„(M(B)). In this case, in view of Eq.(2), the 
phase difference will acquire an extra M(B) dependent contribution and as a result the 
r.h.s. of Eq.(5) will become a nonlinear functional of M(B). The numerical solution of 
this implicit equation for the resulting magnetization ruf = M(B, 5(B))/Mo(0) is shown 
in FigJH for three values of zero-field deficiency parameter 5(0). As is clearly seen, rrif 
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exhibits a field-induced "fishtail" -like behavior typical for underdoped crystals with in- 
tragrain granularity. The extra extremum of the magnetization appears when the applied 
magnetic field B matches an intrinsic chemomagnetic field B^((5(B)) (which now also de- 
pends on B via the above-discussed magnetoconcentration effect). Notice that a "fishtail" 
structure of manifests itself even at zero values of field-free deficiency parameter 5(0) 
(solid line in Fig. 3) thus confirming a field- induced nature of intrinsic nanogranularity 
(Lang et al. 2002, Daeumling et al. 1990, Yang et al. 1993, Gurevich and Pashitskn 1997, 
Moeckley et al. 1993). At the same time, even a rather small deviation from the zero-field 
stoichiometry (with S{0) = 0.1) immediately brings about a paramagnetic Mcissner effect 
at low magnetic fields. Thus, the present model predicts appearance of two interrelated 
phenomena, Meissner paramagnetism at low fields and "fishtail" anomaly at high fields. It 
would be very interesting to verify these predictions experimentally in non-stoichiometric 
superconductors with pronounced networks of planar defects. 

3. MAGNETIC FIELD INDUCED POLARIZATION EFFECTS IN 2D JJA 



In this Section, within the same model of JJAs created by a regular 2D network of 
twin-boundary (TB) dislocations with strain fields acting as an insulating barrier between 
hole-rich domains in underdoped crystals, we discuss charge-related effects which are ac- 
tually dual to the abovc-dcscribcd phase-related chemomagnetic effects. Specifically, wc 
consider a possible existence of a non-zero electric polarization P(5, B) (chemomagnc- 
toelectric effect) and the related change of the charge balance in intrinsically granular 
non-stoichiometric material under the infiuence of an applied magnetic field. In par- 
ticular, we predict an anomalous low-field magnetic behavior of the effective junction 
charge Q(5, B) and concomitant magnetocapacitance C(5, B) in paramagnetic Meissner 
phase and a charge analog of "fishtail" -like anomaly at high magnetic fields along with 
field-tuned weakening of the chemically-induced Coulomb blockade (Sergeenkov 2007). 

Recall that a conventional (zero-field) pair polarization operator within the model 
under discussion reads (Sergeenkov 1997, 2002, 2007) 

JV 

p = ^ g^Xi (6) 

i=l 

In view of Eqs.(l), (2) and (6), and taking into account a usual "phase-number" commu- 
tation relation, [0i,nj] = iSij, it can be shown that the evolution of the pair polarization 
operator is determined via the equation of motion 

dp 1 r o /I 2e T • , / \ /™\ 
^ = [P> = Jij sm 0ij(i)xij- (7) 

Resolving the above equation, we arrive at the following net value of the magnetic-field 
induced longitudinal (along x-axis) electric polarization P{S, B) and the corresponding 
effective junction charge 

T t 

Q(S,B) = ^JdtJdt' J ^sin(/.(x,t')a;e-W/^ (8) 
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where S = 2TTd? is properly defined normalization area, r is a characteristic time (see 
below), and we made a usual substitution j^Ylij ^ ^ I d'^xA{yi,t) valid in the 

long- wavelength approximation (Sergeenkov 2002). 

To capture the very essence of the superconducting analog of the chemomagnetoelectric 
effect, in what follows we assume for simplicity that a stoichiometric sample (with 5 ~ 0) 
does not possess any spontaneous polarization at zero magnetic field, that is P(0, 0) = 0. 
According to Eq.(8), this condition implies = 2nm for the initial phase difference with 
m = 0,±1,±2,... 

Taking the applied magnetic field along the c-axis (and normal to the CuO plane), we 
obtain finally 

Q(^,B) = Qo(^) ^^ + Ml-b-) 

^ ^ ^ ^(l + b2)(l + b2)2 ^ ^ 

for the magnetic field behavior of the effective junction charge in chemically induced 
granular superconductors. Here Qo(^) = eTJo{S)/h with Jo{6) defined earlier, b = B/Bq, 
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Fig. 5: A variation of effective junc- 
tion charge with an applied mag- 
netic field (chemomagnetoelectric ef- 
fect) for different values of oxygen 
deficiency parameter: 5 ~ (solid 
line), 6 = 0.1 (dashed line), and 
S = 0.2 (dotted line). 



b = b — b^, and b^ = B^/Bq — {kBTT/Ti)5 where B^(5) = (/i^r/^)Bo is the chemically- 
induced contribution (which disappears in optimally doped systems with 5 ~ 0), and 
Bq = ^o/wd is a characteristic Josephson field. 

FigOshows changes of the initial (stoichiometric) effective junction charge AQ((5, B) = 
Q(5, B) — Q(5, 0) (solid line) with oxygen deficiency 6. According to Eq.(9), the effective 
charge Q changes its sign at low magnetic fields (driven by non-zero values of S) as soon 
as the chemomagnetic contribution B^(5) exceeds an applied magnetic field B. This is 
nothing else but a charge analog of chemically induced PME. At the same time, Figl6] 
presents a variation of the chemomagnetoelectric effect with concentration (deficiency) for 
different values of the applied magnetic field. Notice that a zero-field contribution (which 
is a true chemoelectric effect) exhibits a maximum around 6c — 0.2, in agreement with 
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Fig. 6: A variation of the chemo- 
magnetoelectric effect with concen- 
tration (deficiency) for different val- 
ues of the apphed magnetic field. 



the classical percolative behavior observed in non-stoichiometric YBa2CusO'j-s samples 
(Gantmakher et al. 1990). 

It is of interest also to consider the magnetic field behavior of the concomitant effective 
flux capacitance C = TdQ,{6, B)/(i$ which in view of Eq.(9) reads 

^/r^N ^ 1 -Sbb-Sb^ + bb^ , , 

C(5,B =Co5 , 10 

(l + b2)(l + b2)3 ^ ^ 

where $ = SB, and Co(5) = rQo(5)/<l>o. 

FigH depicts the behavior of the effective flux capacitance AC(5, B) = 0(5, B) — 
C{6, 0) in applied magnetic field for different values of oxygen deficiency parameter: 5^0 
(solid line), 6 = 0.1 (dashed line), and 6 = 0.2 (dotted line). Notice a decrease of magne- 
tocapacitance amplitude and its peak shifting with increase of S and sign change at low 
magnetic fields which is another manifestation of the charge analog of chemically induced 
PME (Cf. Fig|5]). Up to now, we neglected a possible field dependence of the chemical 
potential /i„ of oxygen vacancies. Recall, however, that in high enough applied magnetic 
fields B, the field-induced change of the chemical potential Afiy(B) = /it,(B) — /i^,(0) be- 
comes tangible and should be taken into account (Abrikosov 1988, Sergeenkov and Ausloos 
1999). As a result, we end up with a superconducting analog of the so-called magneto- 
concentration effect (Sergeenkov 2003) with field induced creation of oxygen vacancies 
Cv(B) = Cv(0) exp(— A/i„(B)//csT) which in turn brings about a "fishtail" -like behavior 
of the high-field chemomagnetization (see Section 2 for more details). FigJHlshows the field 
behavior of the effective junction charge in the presence of the above-mentioned magne- 
toconcentration effect. As it is clearly seen, Q((5(B),B) exhibits a "fishtail" -like anomaly 
typical for previously discussed (Sergeenkov 2003) chemomagnetization in underdoped 
crystals with intragrain granularity. This more complex structure of the effective charge 
appears when the applied magnetic field B matches an intrinsic chemomagnetic field 
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B/B„ 



Fig. 7: The effective ffux capaci- 
tance as a function of applied mag- 
netic field for different values of oxy- 
gen deficiency parameter: 5 ~ 
(solid line), 6 = 0.1 (dashed line), 
and 5 = 0.2 (dotted line). 




B/B„ 



Fig. 8: A "fishtail" -like behavior of 
an effective charge in applied mag- 
netic field in the presence of mag- 
netoconcentration effect (with field- 
induced oxygen vacancies 5(B)) for 
three values of field-free deficiency 
parameter (from top to bottom): 
5(0) ~ (solid line), 5(0) = 0.1 
(dashed line), and 5(0) = 0.2 (dot- 
ted line). 



B^(5(B)) (which now also depends on B via the magnetoconcentration effect). Notice 
that a "fishtail" structure of Q(5(B),B) manifests itself even at zero values of field-free 
deficiency parameter 5(0) (solid line in FiglH]) thus confirming a field-induced nature of 
intrinsic granularity. Likewise, FiglHl depicts the evolution of the effective flux capacitance 
AC(5(B),B) = C(5(B),B) - C(5(0),0) in applied magnetic field B in the presence of 
magnetoconcentration effect (Cf. FigjT]). 
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8(0)=0 
8(0)=0.1 
5(0)=0.2 



Fig. 9: The behavior of the effec- 
tive flux capacitance in apphed mag- 
netic field in the presence of magne- 
toconcentration effect for three val- 
ues of field-free deficiency parame- 
ter: 5{0) ~ (solid line), 5(0) = 0.1 
(dashed line), and 6{0) = 0.2 (dot- 
ted line) . 
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Thus, the present model predicts appearance of two interrelated phenomena dual to 
the previously discussed behavior of chemomagnetizm (see Section 2), namely a charge 
analog of Meissner paramagnetism at low fields and a charge analog of "fishtail" anomaly 
at high fields. To see whether these effects can be actually observed in a real material, 
let us estimate an order of magnitude of the main model parameters. 

Using typical for HTS single crystals values of Al(0) ~ 150nm, d ~ lOnm, and 
jc ~ lO^^A/m?, we arrive at the following estimates of the characteristic Bq ~ 0.5T and 
chemomagnetic B^((5) ~ 0.5Bq fields, respectively. So, the predicted charge analog of 
PME should be observable for applied magnetic fields B < 0.25T. Notice that, for the 
above set of parameters, the Josephson length is of the order of Aj ~ l/^m, which means 
that the small-junction approximation assumed in this paper is valid and the "self- field" 
effects can be safely neglected. 

Furthermore, the characteristic frequencies a; ^ r^^ needed to probe the effects sug- 
gested here are related to the processes governed by tunneling relaxation times r ~ 
h/Jo{d). Since for oxygen deficiency parameter 6 = 0.1 the chemically-induced zero- 
temperature Josephson energy in non-stoichiometric YBCO single crystals is of the order 
of Jo{S) ~ ksTcS ~ ImeV, we arrive at the required frequencies of a; ~ lO^^Hz and at 
the following estimates of the effective junction charge Qo — e = 1.6 x 10~^^C and fiux 
capacitance Cq — lO^^^F. Notice that the above estimates fall into the range of param- 
eters used in typical experiments for studying the single-electron tunneling effects both 
in JJs and JJAs (Makhlin et al. 2001, van Bentum et al. 1988) suggesting thus quite an 
optimistic possibility to observe the above-predicted field induced effects experimentally 
in non-stoichiometric superconductors with pronounced networks of planar defects or in 
artificially prepared JJAs. It is worth mentioning that a somewhat similar behavior of 
the magnetic field induced charge and related fiux capacitance has been observed in 2D 
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electron systems (Chen et al. 1994). 

And finally, it can be easily verified that, in view of Eqs.(6)-(8), the field-induced 
Coulomb energy of the oxygen-depleted region within our model is given by 



with Q(5, B) and C(5, B) defined by Eqs. (9) and (10), respectively. 

A thorough analysis of the above expression reveals that in the PME state (when 
B <^ B^) the chemically- induced granular superconductor is in the so-called Coulomb 
blockade regime (with Ec > Jo), while in the "fishtail" state (for B > B^) the energy 
balance tips in favor of tunneling (with Ec < Jo)- In particular, we obtain that Ec{S, B = 
O.IB^) = |Jo(^) S'-nd Ec{S,B = B^) = |Jo(5). It would be also interesting to check this 
phenomenon of field-induced weakening of the Coulomb blockade experimentally. 

4. GIANT ENHANCEMENT OF THERMAL CONDUCTIVITY IN 2D JJA 

In this Section, using a 2D model of inductive Josephson junction arrays (created 
by a network of twin boundary dislocations with strain fields acting as an insulating 
barrier between hole- rich domains in underdoped crystals), we study the temperature, 
T, and chemical pressure, V/x, dependence of the thermal conductivity (TC) k of an 
intrinsically nanogranular superconductor. Two major effects affecting the behavior 
of TC under chemical pressure arc predicted: decrease of the linear (i.e., VT - in- 
dependent) TC, and giant enhancement of the nonlinear (i.e., VT - dependent) TC 
with [k{T, VT, V/x) - k{T, VT, 0)]/k{T, VT, 0) reaching 500% when chemoelectric field 

— ^V// matches thermoelectric field = S'tVT. The conditions under which these 
effects can be experimentally measured in non-stoichiometric high-Tc superconductors 
are discussed. 

There are several approaches for studying the thermal response of J Js and J JAs based 
on phenomenology of the Josephson effect in the presence of thermal gradients (see, e.g., 
van Harhngen et al. 1980, Guttman et al. 1997, Deppe and Feldman 1994, Sergeenkov 
2002, Sergeenkov 2007 and further references therein). To adequately describe transport 
properties of the above-described chemically induced nanogranular superconductor for all 
temperatures and under a simultaneous infiuence of intrinsic chemical pressure V/i(x) = 
iirf2oVe(x) and applied thermal gradient VT, we employ a model of 2D overdamped 
Josephson junction array which is based on the following total Hamiltonian (Sergeenkov 




(11) 



2002) 



n{t) = nT{t) + nL{t) + n^{t), 



(12) 



where 



N 




(13) 



is the well-known tunneling Hamiltonian, 




(14) 
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accounts for a mutual inductance Lij between grains (and controls the normal state value 
of the thermal conductivity, see below) with = {fi/2e)(j)ij{t) being the total magnetic 

flux through an array, and finally 

N 

n^{t)^J2ni{t)5i^i (15) 

i=l 

describes chemical potential induced contribution with = XjV/x, and being the pair 
number operator. 

According to the above-mentioned scenario, the tunneling Hamiltonian Tirit) in- 
troduces a short-range (nearest-neighbor) interaction between N junctions (which are 
formed around oxygen-rich superconducting areas with phases 0i(t)), arranged in a two- 
dimensional (2D) lattice with coordinates Xi = {xi, Ui). The areas are separated by oxygen- 
poor insulating boundaries (created by TB strain fields e(xij)) producing a short-range 
Josephson coupling J^j = jQ{5)e~^^^^^^'^. Thus, typically for granular superconductors, 
the Josephson energy of the array varies exponentially with the distance Xj^- = Xj — x^ 
between neighboring junctions (with d being an average grain size). The temperature 
dependence of chemically induced Josephson coupling is governed by the following ex- 
pression, Jij{T) = Jij{0)F{T) where 



F(T) = =Wtanh 
^ ' AO 



A(T) ' 
2kBT 



(16) 



and Jij(O) = [A[0) / 2]{Rq / Rij) with A(T) being the temperature dependent gap param- 
eter, Rq = h/Ae^ is the quantum resistance, and Rij is the resistance between grains in 
their normal state. 

By analogy with a constant electric field E, a thermal gradient VT applied to a 
chemically induced JJA will cause a time evolution of the phase difference across insulating 
barriers as follows (Sergeenkov 2002) 

0,,-(i)=(/.°. + a;,,(V//,VT)i (17) 

Here 0^- is the initial phase difference (see below), and a;,,- = 2e(E„ — ET)xj,//i where 
E^ = ^^/^ ^^'^ = 5'rVT are the induced chemoelectric and thermoelectric fields, 
respectively. St is the so-called thermophase coefficient (Sergeenkov 1998a) which is 
related to the Seebeck coefficient Sq as follows, St = {l/d)So (where / is a relevant 
sample's size responsible for the applied thermal gradient, that is |VT| = AT//). 

We start our consideration by discussing the temperature behavior of the conventional 
(that is linear) thermal conductivity of a chemically induced nanogranular superconductor 
paying a special attention to its evolution with a mutual inductance Lij. For simplicity, in 
what follows we limit our consideration to the longitudinal component of the total thermal 
flux Q(t) which is defined (in a q-spacc representation) via the total energy conservation 
law as follows 



Q{t) = lim 

q^O 



where — dHd/ dt with 



(18) 

n^{t) = d'^xe''^''n{x,t) (19) 
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Here s = 27rci^ is properly defined normalization area, and we made a usual substitution 
Tf ~^ 1 1 (PxA{yi, t) valid in the long-wavelength approximation (q ^ 0). 

In turn, the heat flux Q,{t) is related to the Zmear thermal conductivity (LTC) tensor 
i'^ai3 by the Fourier law as follows (hereafter, {a, /3} = x, y, z) 



d< Qa > 

a(V^T) 



(20) 

VT=0 



where 

< 



Qa > = - / dt< Qa{t) > (21) 
T Jo 

Here V is sample's volume, r is a characteristic Josephson tunneling time for the network, 
and < ... > denotes the thermodynamic averaging over the initial phase differences 0° 

< M<I>1) >= I n d<l>lA{<f>l)e-^^o (22) 



with an effective Hamiltonian 



Horn - f (23) 

Jo ^ J ^ 

Here, jS = l/Zc^T, and Z = J^^ Ylij d(j)^^-e~^^'^ is the partition function. The above-defined 
averaging procedure allows us to study the temperature evolution of the system. 

Taking into account that in JJAs (Eichenberger et al. 1996) Lij oc Rij, we obtain 
Lij — Loexp(| Xjj \/d) for the explicit x-dependence of the weak- link inductance in our 
model. Finally, in view of Eqs.(12)-(23), and making use of the usual "phase- number" 
commutation relation, [(j)i,nj\ = iSij, wc find the following analytical expression for the 
temperature and chemical gradient dependence of the electronic contribution to linear 
thermal conductivity of a granular superconductor 

Ka(3{T, V//) = >^o[SaMT^: e) + /?L(T)i/(T, e)/.^^)] (24) 

where 

fap{e) = I [SapA{e) - e^epB{e)] (25) 

with 

^(^)=(^^ + ;tan-^ (26) 

and 

,2(1^,2)3 +^tan-^^ (27) 

Here, kq = Nd'^Sr^/ VLo, /^^(T) = 27r/c(T)Lo/$o with /c;(T) = {2e/h)J{T) being the 
critical current; e = + + el with = E^/Eq where Eq = h/2edT is a characteristic 
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field. In turn, the above-introduced "order parameters" of tlie system, ri{T,e) 



and z/fT, e) =< sin 



6?- >, are defined as follows 



71 4 



TT— (2n+l)2 



and 



where 



/5J(T) 



sinh 13^ 



1 1 



H — tan ^ e 



(28) 

(29) 
(30) 



Here J(T) is given by Eq.(17), and In{x) stand for the modified Bessel functions. 

Turning to the discussion of the obtained results, we start with a more simple zero- 
pressure case. The relevant parameters affecting the behavior of the LTC in this particular 
case include the mutual inductance Lq and the normal state resistance between grains Rn- 
For the temperature dependence of the Josephson energy (see Eq.(17)), we used the well- 
known (Sergeenkov 2002) approximation for the BCS gap parameter, valid for all temper- 



atures, A(T) = A(0) tanh 7 



with 7 = 2.2. Despite a rather simplified nature of 
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Fig. 10: Temperature dependence of the zero-pressure (Vfi = 0) linear thermal conduc- 
tivity for r„ = 0.1 (left) and r„ = 1 (right) for different values of the SQUID parameter 
(from bottom to top): /3l(0) = 1, 3, 5, 7, and 9. 

our model, it seems to quite reasonably describe the behavior of the LTC for all tempera- 
tures. Indeed, in the absence of intrinsic chemical pressure (V/x = 0), the LTC is isotropic 
(as expected), Ka/3(T,0) = 5„/3Kl(T,0) where kl{T,0) = ko[??(T, 0) + 2/3l(T)i/(T, 0)] van- 
ishes at zero temperature and reaches a normal state value = ^^.(Tc, 0) = (7r/2)Ko at 
T = Tc. FigJTOl shows the temperature dependence of the normalized LTC ^^(T, 0)/k„ 
for different values of the so-called SQUID parameter /3l(0) = 27r/c(0)Lo/$o (increasing 
from the bottom to the top) and for two values of the resistance ratio r„ = Rq/ = 0.1 
and Tn = Ro/Rn = 1. First of all, with increasing of the SQUID parameter, the LTC 
evolves from a fiat-like pattern (for a relatively small values of Lq) to a low-temperature 
maximum (for higher values of (3l{0)). Notice that the peak temperature Tp is practi- 
cally insensitive to the variation of inductance parameter Lq while being at the same time 
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strongly influenced by resistivity Rn- Indeed, as it is clearly seen in FigJTO| a different 
choice of r„ leads to quite a tangible shifting of the maximum. Namely, the smaller is 
the normal resistance between grains i?„ (or the better is the quality of the sample) the 
higher is the temperature at which the peak is developed. As a matter of fact, the peak 
temperature Tp is related to the so-called phase-locking temperature Tj (which marks 
the establishment of phase coherence between the adjacent grains in the array and always 
lies below a single grain superconducting temperature Tc) which is usually defined via 
an average (per grain) Josephson coupling energy as J(Tj,r„) = /c^Tj. Indeed, it can 
be shown analytically that for Tj < T < Tc, Tj(r„) ~ r„Tc. 

Turning to the discussion of the LTC behavior under chemical pressure, let us assume, 
for simplicity, that V/i = (V^/x, 0, 0) with oxygen-deficiency parameter 5 controlled chem- 
ical pressure V^^/i ~ ^v{S)/d, and VT = (V^T, Vj^T, 0). Such a choice of the external 
fields allows us to consider both parallel K2;2;(T, V/i) and perpendicular Kj^y(T, V/i) com- 
ponents of the LTC corresponding to the two most interesting configurations, V/i||VT 
and V/i-LVT, respectively. FigHT] demonstrates the predicted chemical pressure depen- 
dence of the normalized LTC Akl(T, V/i) = ^^.(T, V/i) — kl(T, 0) for both configurations 
taken at T = 0.2Tc (with r„ = 0.1 and /9l(0) = 1). First of all, we note that both com- 
ponents of the LTC are decreasing with increasing of the pressure E^/Eq = fiv{S)T/h. 
And secondly, the normal component Kyy decreases more slowly than the parallel one 
Kxx, suggesting thus some kind of anisotropy in the system. In view of the structure 
of Eq.(25), the same behavior is also expected for the temperature dependence of the 
chemically-induced LTC, that is AKi(T, V/i)/Ki(T, 0) < for all gradients and tem- 
peratures. In terms of the absolute values, for T = 0.2Tc and E^ = Eg, we obtain 
[AKi(T, V/i)/KL(T,0)]^^ = 90% and [Akl(T, V/i)/KL(T, 0)]yj, = 60% for attenuation of 
LTC under chemical pressure. Let us turn now to the most intriguing part of this Section 




Fig. 11: The dependence of the 
linear thermal conductivity on 
the chemical pressure for paral- 



^ E /E ^ lei (V/i||VT) and perpendicular 

(V/i-LVT) configurations. 

and consider a nonlinear generalization of the Fourier law and very unusual behavior of 
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the resulting nonlinear thermal conductivity (NLTC) under the influence of chemical pres- 
sure. In what follows, by the NLTC we understand a VT-dependent thermal conductivity 
rejf (T, V//) = Kap{T, V//; VT) which is defined as follows 



«:2f(T,V/.)^-- 



d< Qa > 

a(V^T) 



(31) 



VT^O 



with < Qa > given by Eq.(21). 

Repeating the same procedure as before, we obtain finally for the relevant components 
of the NLTC tensor 



where 
with 

and 



</r(T, V/x) = Ko[5apv{T, eeff) + /3l(T)i/(T, eeff)D^pieeff)], 



, 3 + lle^ - lle^ - 3e^ 3 
C{e) — — — tan e 



(32) 

(33) 

(34) 

(35) 



e4(l + e2)4 

Here, e^^^ = e° - e?; where e° = E^^/Eq and = E^/Eo with E^ = 5tV«T; other 
parameters (77, v, B and fap) are the same as before but with e — > ee//. As expected, in 
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Fig. 12: The dependence of the non- 
linear thermal conductivity on the 
chemical pressure for different values 
of the applied thermal gradient ey = 
StVT/Eo {er = 0.2,0.4,0.6,0.8, and 
1.0, increasing from bottom to top). 



the limit (or when E^ » Er), from Eq.(32) we recover all the results obtained 

in the previous section for the LTC. Let us see now what happens when thermoelectric 



19 



field Et = StVT becomes comparable with chemoelectric field E^. Fig{T2] depicts the 
resulting chemical pressure dependence of the parallel component of the NLTC tensor 
Ak^^^(T,E^) = Kf^^(T,E^) - Kf/'(T,0) for different values of the dimensionless pa- 
rameter eT = E7-/E0 (the other parameters are the same as before). As is clearly seen 
from this picture, in a sharp contrast with the pressure behavior of the previously con- 
sidered LTC, its nonlinear analog evolves with the chemoelectric field quite differently. 
Namely, NLTC strongly increases for small pressure values (with E^ < E^,), reaches a 
pronounced maximum at E^ = E^ = f Et, and eventually declines at higher values of ^ 
(with E^ > Em). Furthermore, as it directly follows from the very structure of Eq.(32), 
a similar "reentrant-like" behavior of the non/mear thermal conductivity is expected for 
its temperature dependence as well. Even more remarkable is the absolute value of the 
pressure- induced enhancement. According to FigJT2| it is easy to estimate that near 
maximum (with E^ = E^ and Et = Eq) one gets Ak^^^{T, E^)/Kf/'(T, 0) ~ 500%. 

To understand the above-obtained rather unusual results, let us take a closer look 
at the chemoelectric field induced behavior of the Josephson voltage in our system (see 
Eq.(17)). Clearly, strong heat conduction requires establishment of a quasi- stationary 
(that is nearly zero-voltage) regime within the array. In other words, the maximum of 
the thermal conductivity under chemical pressure should correlate with a minimum of 
the total voltage in the system, V{Vfi) = (^) < >= Vo{e - e^) where e = E^/Eq 

and Vo = Eqc? = h/2eT is a characteristic voltage. For linear TC (which is valid only for 
small thermal gradients with st = Et/Eq -C 1), the average voltage through an array 
Vl{V^) ~ Vo(E^/Eo) has a minimum at zero chemoelectric field (where LTC indeed has 
its maximum value, see Fig. 11) while for nonlinear TC (with — 1) we have to consider 
the total voltage V^(V/i) which becomes minimal at E^ = E^ (in a good agreement with 
the predictions for NLTC maximum which appears at E^ = fE^, see Fig JT^ . 

To complete our study, let us estimate an order of magnitude of the main model 
parameters. Starting with chemoelectric fields E^ needed to observe the above-predicted 
nonlinear field effects in nanogranular superconductors, we notice that according to FigJT2l 
the most interesting behavior of NLTC takes place for E^ ~ Eq. Using typical YBCO 
parameters, eo = 0.01, Qq = with = 0.2nm, and K = 115GPa, we have /j,^ = 
eoKQo ~ ImeV for an estimate of the chemical potential in HTS crystals, which defines 
the characteristic Josephson tunneling time r ~ ^//i^, ~ 5 x 10~^^s and, at the same 
time, leads to creation of excess vacancies with concentration Cy = e"'^"/'^'^'^ ~ 0.75 at 
T = 0.2Tc (equivalent to a deficiency value of 5 ~ 0.25). Notice that in comparison 
with this linear defects mediated channeling (osmotic) mechanism, the probability of the 
conventional oxygen diffusion in these materials D oc e~^'^/^'^'^ (governed by a rather high 
activation energy Ud — leV) is extremely low under the same conditions (D -C 1). 

Furthermore, taking d ^ lOnm for typical values of the average "grain" size (cre- 
ated by oxygen-rich superconducting regions), we get Eq = h/2edT ~ 5 x l{f'V/m and 
|V/i| = ^v/d ~ 10^e\^/m for the estimates of the characteristic field and chemical poten- 
tial gradient (intrinsic chemical pressure), respectively. On the other hand, the maximum 
of NLTC occurs when this field nearly perfectly matches an "intrinsic" thermoelectric field 
Ey = StVT induced by an applied thermal gradient, that is when E^ ^ Eq — Ey. Recall- 
ing that St = {l/d)So and using So — 0.5^V/K and / ^ 0.5mm for an estimate of the lin- 
ear Seebeck coefficient and a typical sample's size, we obtain VT ~ Eo/S*]" — 2 x 10^-ft"/m 
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for the characteristic value of apphed thermal gradient needed to observe the predicted 
here giant chemical pressure induced effects. Let us estimate now the absolute value of 
the linear thermal conductivity governed by the intrinsic Josephson junctions. Recall 
that within our model the scattering of normal electrons is due to the presence of mutual 
inductance between the adjacent grains Lq which is of the order of Lq ~ fi^d ~ IfH as- 
suming d = lOnm for an average "grain" size. In the absence of chemical pressure effects, 
the temperature evolution of LTC is given by kl(T,0) = Ko[r]{T,0) + 2/5l(T)z/(T, 0)] 
where kq = Nd'^Sx^o/VLQ. Assuming V ~ Nd"^! for the sample's volume, using the 
above-mentioned expression for St, and taking /?l(0) = 5 and r„ = 0.1 for the value of 
the SQUID parameter and the resistance ratio, we obtain k,l{0.2Tc, 0) — IW/mK for an 
estimate of the maximum of the LTC (see Fig JTOl) . 

And finally, it is worth comparing the above estimates for inductively coupled grains 
(Sergeenkov 2002) with the estimates for capacitively coupled grains (Sergeenkov 2007) 
where the scattering of normal electrons is governed by the Stewart-McCumber parameter 
/5(7(T) = 27rJc(T)Co-R,^/$o due to the presence of the normal resistance i?„ and mutual 
capacitance Cq between the adjacent grains. The latter is estimated to be Cq ~ laF 
using d = lOnm for an average "grain" size. Furthermore, the critical current /c(0) can 
be estimated via the critical temperature Tc as follows, /c(0) — 27ikBTc/^o which gives 
/c(0) ~ lOfxA (for Tc ~ 90K) and leads to PdO) ^ 3 for the value of the Stewart- 
McCumber parameter assuming i?„ ~ Rq for the normal resistance which, in turn, results 
in g ~ $o/-Rn — 10"^^C and Ec = /2Cq ~ Q.leV for the estimates of the "grain" charge 
and the Coulomb energy. Using the above-mentioned expressions for Sq and /3c (0), we 
obtain ~ lO'^W/mK for the maximum of the capacitance controlled LTC which is 
actually much smaller than a similar estimate obtained above for inductance controlled 
kl (Sergeenkov 2002) but at the same time much higher than phonon dominated heat 
transport in granular systems (Deppe and Feldman 1994). 

5. THERMAL EXPANSION OF A SINGLE JOSEPHSON CONTACT AND 2D JJA 



In this Section, by introducing a concept of thermal expansion (TE) of a Josephson 
junction as an elastic response to an effective stress field, we study (both analytically 
and numerically) the temperature and magnetic field dependence of TE coefficient a in 
a single small junction and in a square array. In particular, we found (Sergeenkov et 
al. 2007) that in addition to field oscillations due to Fraunhofer-like dependence of the 
critical current, a of a small single junction also exhibits strong flux driven temperature 
oscillations near Tc. We also numerically simulated stress induced response of a closed 
loop with finite self-inductance (a prototype of an array) and found that a of a 5 x 5 
array may still exhibit temperature oscillations if the applied magnetic field H is strong 
enough to compensate for the screening induced effects. 

Since thermal expansion coefficient a(T, H) is usually measured using mechanical 
dilatometers (Nagel et al. 2000), it is natural to introduce TE as an elastic response of 
the Josephson contact to an effective stress field a (D'yachenko et al. 1995, Sergeenkov 
1998b, Sergeenkov 1999). Namely, we define the TE coefficient (TEC) a(T, H) as follows: 

«(T,H) = ^ (36) 
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where an appropriate strain field e in the contact area is related to the Josephson energy 
Ej as follows {V is the volume of the sample): 



1 



dEj 
da 



(37) 

(T=0 



For simplicity and to avoid self-field effects, we start with a small Josephson contact 
of length w < Xj (Aj — y^$o/ l^odjc is the Josephson penetration depth) placed in a 
strong enough magnetic field (which is applied normally to the contact area) such that 
H > ^0/211 Xjd, where d — 2XL + t, Xl is the London penetration depth, and t is an 
insulator thickness. 

The Josephson energy of such a contact in applied magnetic field is governed by a 
Praunhofer-hke dependence of the critical current (Orlando and Delin 1991): 

Ej^j(^l-^-^cosifoy (38) 

where </? = 7r$/$o is the frustration parameter with $ = Hwd being the flux through the 
contact area, ipo is the initial phase difference through the contact, and J oc e~*/^ is the 
zero-field tunneling Josephson energy with ^ being a characteristic (decaying) length and 
t the thickness of the insulating layer. The self-field effects (screening), neglected here, 
will be considered later for an array. 

Notice that in non-zero applied magnetic field H, there are two stress-induced contri- 
butions to the Josephson energy Ej, both related to decreasing of the insulator thickness 
under pressure. Indeed, according to the experimental data (D'yachenko ct al. 1995), the 
tunneling dominated critical current Ic in granular high-Tc superconductors was found 
to exponentially increase under compressive stress, viz. /c(c) = /c(0)e'*'^. More specif- 
ically, the critical current at cr = 9kbar was found to be three times higher its value at 
(7 = l.Bkbar, clearly indicating a weak-links-mediated origin of the phenomenon. Hence, 
for small enough a we can safely assume that (Sergeenkov 1999) t(a) ~ t(0)(l — Pa/ao) 
with (To being some characteristic value (the parameter f3 is related to the so-called ul- 
timate stress am SiS (3 = ao/am)- As a result, we have the following two stress- induced 
effects in Josephson contacts: 

(I) amplitude modulation leading to the explicit stress dependence of the zero-field 
energy 

J(T, a) = J(T, 0)6^'"/'"° (39) 

with 7 = (3t{0)/^, and 

(II) phase modulation leading to the explicit stress dependence of the fiux 

$(T, H, a) = Hw(i(T, a) (40) 

with 

d{T, a) = 2Al(T) + t(0)(l - Pa/ao) (41) 

Finally, in view of Eqs.(36)-(41), the temperature and field dependence of the small 
single junction TEC reads (the initial phase difference is conveniently fixed at ipo = n): 

a(T, H) = a(T, 0) [1 + F(T, H)] + e(T, 0) ^^^^^^^ (42) 
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where 



F(T,H) 



sill i 



C / sin V5 



with 



7r$(T,H,0) H 



COS(^ 



(43) 

(44) 

(45) 
and 

Here, Ho(T) = ^o/nwd{T, 0) with d{T, 0) = 2Al(T) +i(0). For the exphcit temperature 



a{T,0) 



$0 Ho(T) 
_ c^e(T,0) 
~ dT 



o 

to 



CO 
< 




T/T, 



Fig. 13: Temperature dependence of 
the flux driven strain field in a single 
short contact for different values of 
''■^ the frustration parameter f accord- 
ing to Eqs.(36)-(48). 



dependence of J(T, 0) = $o-^c(T)/27r we use the well-known (Meservey and Schwartz 
1969, Sergeenkov 2002) analytical approximation of the BCS gap parameter (valid for all 

temperatures), A(T) = A(0) tanh ^2.2^^^^ with A(0) = 1.76A;bTc which governs 

the temperature dependence of the Josephson critical current 



7c(T) = lc{0) 



A(T) 
A(0) 



tanh 



' A(T) ' 
2kBT 



(47) 



while the temperature dependence of the London penetration depth is governed by the 
two-fiuid model: 

Al(T) = , (48) 
' VI - (T/Tc)2 ^ ' 

From the very structure of Eqs.(36)-(44) it is obvious that TEC of a single contact will 
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Fig. 14: Temperature dependence of 
flux driven normalized TEC in a sin- 
gle small contact for different val- 
ues of the frustration parameter f 
(for the same set of parameters as in 
Fig. 13) according to Eqs.(36)-(48). 



exhibit /le/d oscillations imposed by the Fraunhofer dependence of the critical current Ic- 
Much less obvious is its temperature dependence. Indeed, Fig{T3] presents the temperature 
behavior of the contact area strain field Ae(T,f) = e(T,f) - e(T, 0) (with t{0)/^ = 
1, ^/Xl{0) = 0.02 and (3 = 0.1) for different values of the frustration parameter f = 
H/Ho(0). Notice characteristic flux driven temperature oscillations near Tc which are 
better seen on a semi-log plot shown in Fig{T3] which depicts the dependence of the 
properly normalized field-induced TEC A«(T,f) = a(T,f) — a(T,0) as a function of 
1 — T/Tc for the same set of parameters. 

To answer an important question how the neglected in the previous analysis screening 
effects will affect the above-predicted oscillating behavior of the field-induced TEC, let us 
consider a more realistic situation with a junction embedded into an array (rather than an 
isolated contact) which is realized in artificially prepared arrays using photolithographic 
technique that nowadays allows for controlled manipulations of the junctions parameters 
(Newrock et al. 2000). Besides, this is also a good approximation for a granular super- 
conductor (if we consider it as a network of superconducting islands connected with each 
other via Josephson links). Our goal is to model and simulate the elastic response of such 
systems to an effective stress a. For simplicity, we will consider an array with a regular 
topology and uniform parameters (such approximation already proved useful for describ- 
ing high-quality artificially prepared structures, see, e.g., Sergeenkov and Araujo-Moreira 



Let us consider a planar square array as shown in FigJTSl The total current includes 
the bias current flowing through the vertical junctions and the induced screening currents 
circulating in the plaquette (Nakajima and Sawada 1981). This situation corresponds to 
the inclusion of screening currents only into the nearest neighbors, neglecting thus the 
mutual inductance terms (Phillips et al. 1993). Therefore, the equation for the vertical 



2004). 
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contacts will read (horizontal and vertical junctions are denoted by superscripts h and f , 
respectively) : 




Fig. 15: Left: sketch of a regular square array (a single plaquette). Right: electrical 
scheme of the array with the circulating currents. The bias current is fed via virtual loops 
external to the array. 

+ TT^^ + ^cSinC,, = A/? . + /, (49) 



2e cit2 2ei? dt 



where A//^- = /f^ — and the screening currents V obey the fluxoid conservation 

condition: 

- + - + = - (50) 

Recall that the total flux has two components (an external contribution and the contri- 
bution due to the screening currents in the closed loop) and it is equal to the sum of the 
phase differences describing the array. It is important to underline that the external flux 
in Eq. (!50!) . ?7 = 2'7r<l>'^^*/'^05 is related to the frustration of the whole array, i.e., this is the 
flux across the void of the network (Araujo-Moreira et al. 1997, Araujo-Moreira et al. 
2005, Grimaldi et al. 1996), and it should be distinguished from the previously introduced 
applied magnetic field H across the junction barrier which is related to the frustration of a 
single contact f = 27rH(ity/<l>o and which only modulates the critical current /c(T,H,cr) 
of a single junction while inducing a negligible flux into the void area of the array. 

For simplicity, in what follows we will consider only the elastic effects due to a uniform 
(homogeneous) stress imposed on the array. With regard to the geometry of the array, 
the deformation of the loop is the dominant effect with its radius a deforming as follows: 

a(a) = ao(l - xo"/cro) (51) 

As a result, the self- inductance of the loop L[a) = fiQaF{a) (with F{a) being a geometry 
dependent factor) will change accordingly: 

L{a) = Lo(l - Xg(r/ao) (52) 

The relationship between the coefficients x ^i-nd Xg is given by 

Xg = {1 + aoBg) X (53) 
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0.0 0.3 0.6 0.9 
f 



Fig. 16: Numerical simulation re- 
sults for an array 5x5 (solid line) 
and a small single contact (dashed 
line). The dependence of the nor- 
malized TEC on the frustration pa- 
rameter f (applied magnetic field H 
across the barrier) for the reduced 
temperature T/Tc = 0.95. The pa- 
rameters used for the simulations: 
r] = 0,P = 0.1, m/^ = 1, e/AL = 
0.02, (3l = 10, 76 = 0.95, and Xg = 
X = 0.01. 



where Bg = (4^) . It is also reasonable to assume that in addition to the critical 

a F{a) \ da / ao 

current, the external stress will modify the resistance of the contact: 

as well as capacitance (due to the change in the distance between the superconductors): 

C{a) = ^° ^ Co(l + x^/^o) (55) 

To simplify the treatment of the dynamic equations of the array, it is convenient to 
introduce the standard normalization parameters such as the Josephson frequency: 



27r/c(0) 



the analog of the SQUID parameter: 



and the dissipation parameter: 



& ^ (57) 

<i>0 



2t:Ic{Q)CoRI 




Combining Eqs. (l49l) and (!50l) with the stress- induced effects described by Eqs. (!54l) and 
(1551) and using the normalization parameters given by Eqs. (l56l) - (j58l) . we can rewrite the 
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equations for an array in a rather compact form. Namely, the equations for vertical 
junctions read: 



1 e-x^l'^^ 



,l^^ex<^l-^ sin0^^^. = 7,+ 



/3l (1 - XsC^M) 



Here an overdot denotes the time derivative with respect to the normalized time (inverse 
Josephson frequency), and the bias current is normalized to the critical current without 
stress, 76 = hllc{^). 

The equations for the horizontal junctions will have the same structure safe for the 
explicit bias related terms: 



1 e-^"!'"'' : 



^ + ^ (j>^. + e^"/"° sin (^^ 



Pl (1 - XgCr/c^o) 



Finally, Eqs. (!59i) and ( 1601) should be complemented with the appropriate boundary condi- 
tions (Binder et al. 2000) which will include the normalized contribution of the external 
flux through the plaquette area rj = 27r$'^^'*/$o- It is interesting to notice that Eqs. fl59l) 
and fl60l) will have the same form as their stress-free counterparts if we introduce the 
stress-dependent renormalization of the parameters: 



/3L=/3L(l-xW^o)e^'^/^« 



r] = r]{l- 2xa/ao) 



lb = 7be 



-XO"/o"0 



(62) 
(63) 
(64) 
(65) 



Turning to the discussion of the obtained numerical simulation results, it should be 
stressed that the main problem in dealing with an array is that the total current through 
the junction should be retrieved by solving self-consistently the array equations in the 
presence of screening currents. Recall that the Josephson energy of a single junction for 
an arbitrary current / through the contact reads: 



Ej{T,i,I) = Ej{T,f,Ic) 




(66) 
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Fig. 17: Numerical simulation results for an array 5x5. The influence of the flux across 
the void of the network rj frustrating the whole array on the temperature dependence of 
the normalized TEC for different values of the barrier fleld f frustrating a single junction 
for 7b = 0.5 and the rest of parameters same as in FigJTGl 



The important consequence of Eq. fl66|) is that if no current flows in the array's junction, 
such junction will not contribute to the TEC (simply because a junction disconnected 
from the current generator will not contribute to the energy of the system). 

Below we sketch the main steps of the numerical procedure used to simulate the stress- 
induced effects in the array: 
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(1) a bias point h is selected for the whole array; 

(2) the parameters of the array (screening, Josephson frequency, dissipation, etc) are 
selected and modified according to the intensity of the applied stress a; 



(3) the array equations are simulated to retrieve the static configuration of the phase 
differences for the parameters selected in step 2; 

(4) the total current flowing through the individual junctions is retrieved as: 

i:f = Icsm^f (67) 



(5) the energy dependence upon stress is numerically estimated using the value of the 



total current I^j (which is not necessarily identical for all junctions) found in step 
4 via Eq.(l67D; 

(6) the array energy Ef is obtained by summing up the contributions of all junctions 



with the above-found phase differences 



(7) the stress-modified screening currents //j(T, H, a) are computed using Eq. (!50l) and 
inserted into the magnetic energy of the array E^j = j^Hij^I- jY; 

(8) the resulting strain field and TE coefficient of the array are computed using numer- 
ical derivatives based on the finite differences: 



1 

V 



A {Et, + Ej) 



(68) 



Ao-^O 



a(T,H)^^ (69) 

The numerical simulation results show that the overall behavior of the strain field and 
TE coefficient in the array is qualitatively similar to the behavior of the single contact. 
In Fig{T6] we have simulated the behavior of both the small junction and the array as a 
function of the field across the barrier of the individual junctions in the presence of bias 
and screening currents. As is seen, the dependence of a(T, f) is very weak up to f ~ 0.5, 
showing a strong decrease of about 50% when the frustration approaches f = 1. 

A much more profound change is obtained by varying the temperature for the fixed 
value of applied magnetic field. Fig{T7| depicts the temperature behavior of a(T, f ) (on 
semi-log scale) for different field configurations which include barrier field / frustrating a 
single junction and the flux across the void of the network rj frustrating the whole array. 
First of all, comparing Fig JTTT a) and Fig{T3]we notice that, due to substantial modulation 
of the Josephson critical current /c(T, H) given by Eq. flSSl) . the barrier field f has similar 
effects on the TE coefficient of both the array and the single contact including temperature 
oscillations. However, finite screening effects in the array result in the appearance of 
oscillations at higher values of the frustration f (in comparison with a single contact). On 
the other hand, FigIT7l(b-d) represent the influence of the external field across the void 
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1] on the evolution of a(T,f). As is seen, in comparison with a field-free configuration 
(shown in FigJTWa)). the presence of external field r] substantially reduces the magnitude 
of the TE coefficient of the array. Besides, with rj increasing, the onset of temperature 
oscillations markedly shifts closer to Tc- 

6. SUMMARY 

In this Chapter, using a realistic model of 2D Josephson junction arrays (created 
by 2D network of twin boundary dislocations with strain fields acting as an insulating 
barrier between hole-rich domains in underdoped crystals), we considered many novel 
effects related to the magnetic, electric, elastic and transport properties of Josephson 
nanocontacts and nanogranular superconductors. Some of the topics covered here include 
such interesting phenomena as chemomagnetism and magnetoelectricity, electric analog 
of the "fishtai" anomaly and field-tuned weakening of the chemically-induced Coulomb 
blockade as well as a giant enhancement of nonlinear thermal conductivity (reaching 500% 
when the intrinsically induced chemoelectric field Efj_ oc |V/i|, created by the gradient of 
the chemical potential due to segregation of hole producing oxygen vacancies, closely 
matches the externally produced thermoelectric field Et oc |VT|). Besides, we have 
investigated the infiuence of a homogeneous mechanical stress on a small single Josephson 
junction and on a plaquette (array of 5 x 5 junctions) and have shown how the stress- 
induced modulation of the parameters describing the junctions (as well as the connecting 
circuits) produces such an interesting phenomenon as a thermal expansion (TE) in a single 
contact and two-dimensional array (plaquette). We also studied the variation of the TE 
coefficient with an external magnetic field and temperature. In particular, near Tc (due 
to some tremendous increase of the effective "sandwich" thickness of the contact) the field- 
induced TE coefficient of a small junction exhibits clear temperature oscillations scaled 
with the number of fiux quanta crossing the contact area. Our numerical simulations 
revealed that these oscillations may actually still survive in an array if the applied field 
is strong enough to compensate for finite screening induced self-field effects. 

The accurate estimates of the model parameters suggest quite an optimistic possibility 
to experimentally realize all of the predicted in this Chapter promising and important for 
applications effects in non-stoichiometric nanogranular superconductors and artificially 
prepared arrays of Josephson nanocontacts. 
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